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Subgap resonances in the I-V curves of layered superconductors are explained by the coupling between Joseph- 
son oscillations and phonons with dispersion in c-direction. In the presence of a magnetic field applied parallel 
to the layers additional structures due to fluxon motion appear. Their coupling with phonons is investigated 
theoretically and a shift of the phonon resonances in strong magnetic fields is predicted. 
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1. Excitation of phonons by Josephson os- 
cillations 

The c-axis transport in the highly anisotropic 
cuprate-superconductors Tl 2 Ba 2 Ca 2 Cu30io+,5 
(TBCCO) and Bi 2 Sr 2 CaCu 2 8+(5 (BSCCO) can 
well be described by a model where the super- 
conducting Cu0 2 layers are coupled by tunneling 
barriers forming a stack of Josephson junctions 

Subgap structures [gJJ] in the I-V characteris- 
tics could be explained by the excitation of longi- 
tudinal c-axis phonons by the Josephson oscilla- 
tions in resistive junctions (gj-Q. In Ref. || a mi- 
croscopic theory for the coupling between Joseph- 
son oscillations and phonons has been derived 
for the case of a short stack of Josephson junc- 
tions where it can be assumed that the gauge- 
invariant phase difference j n (x,t) — ip n (x,t) — 
ip n +i{x, t) — -# f dzA z (x, z, t) between super- 
conducting layers n and n + 1 is constant along 
the layers (independent of x). In the simplest 
model the tunneling current through a junction 
n is relatexd the bias current density jb by the 
RSJ equation 



j b = j c sin7„(t) + aE n (t) + D n (t). 



(1) 



Here E n (t) is the average electric field across 
the junction which is related to the phase differ- 



ence 7„ by the second Josephson equation %-j n = 
2edE n where d is the thickness of the barrier. 

D n = E£ is the field generated by the oscil- 
lating conduction electron charges on the Cu0 2 
layers. It can be expressed by the average field 
in the barrier and the ionic polarization, D n = 
eoE n + P n . The polarization P n is proportional 
to the lattice displacements of ions. Phonons are 
excited by the field of the oscillating electronic 
charges on the Cu0 2 layers. The field D(k z ) can 
be expressed by a generalized dielectric function 
D{k z ) = e e^(k z ,uj)E(k z ) of the form 



e Ph 



(a„ W )=(i-x; 



\n(k z ,\)\ 



uj 2 (k z , A) — lu 2 



(2) 



Here uj(k z ,\) are the eigenfrequencies of the 
dynamical matrix (including long-range Coulomb 
forces). The oscillator strength \il(k z , A)| 2 takes 
care of the fact that ions inside the barrier and on 
the Cu0 2 layers are excited by different fields and 
have to be counted differently in the polarization 
in the RSJ-equation (jl]). 

The dielectric function has zeros at the frequen- 
cies of longitudinal c-axis phonons. It can also be 
written in the more common form 



ph 



{k z ,u) 



v fe,A)| ; 
2^«a(fc aj A)-, 



(3) 



where we have included a background DK. The 



frequencies u>(k z , A) where the dielectric function 
has poles can also be calculated from a dynami- 
cal matrix where long-range Coulomb forces have 
been subtracted. In the limit k z — > they corre- 
spond to the transversal optical eigenfrequencies 
of the system. 

In the resistive state the phase difference of bar- 
rier n has the form 



7»(«) 



■Ut+ 5^ n (t) 



(4) 



where 5j n (t) oscillates with the same Josephson 
frequency u>, while for a barrier in the supercon- 
ducting state the term cut is missing. Inserting 
this ansatz in the RSJ equations we can solve for 
the dc part and the oscillating part. In the special 
case where only one junction of the stack is in the 
resistive state we obtain for the normalized dc- 
current Jb := jb/jc as function of the dc-voltage 
V = hu>/(2e): 



Jb = Jqp(V) + 



jj e 2 + a/(e uj) 
2^el + (e 2 +a/(e Lu)r- 



(5) 



Here e(u>) = ei(u;) + it2{oj) is an averaged phonon 
dielectric function defined by 



e(w)=G- 1 (w) + ^f- 
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r ph 
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(6) 



(7) 



Here o)j = w}^/l — (j/j c ) 2 and w} = 2edj c /(heo) 
is the bare Josephson plasma frequency. In the 
case of a dispersionless phonon we just have 



e(u>) 



— P ph 



(«). 



In the general case the I-V curve shows peaks 
at the zeros of the real part of e(u>) corresponding 
to c-axis phonons with a high density of states 
and non- vanishing oscillator strength |f2| 2 . Be- 
sides optical T-point phonons also phonons from 
the edge of the Brillouin zone at k z = w/d 
and even acoustical phonons may contribute to 
phonon resonances in the I-V curves g. 

The coupling to an acoustical phonon at 
k z — w/d may explain a resonance observed for 
TBCCO at 3.2 mV by Seidel et al. § at a fre- 
quency/voltage which is lower then any optical 
phonon branch expected from model calculations. 



A fit to the experimental data can be made with 
reasonable values of the oscillator strength and 
phonon damping which are compatible with op- 
tical experiments. Also a double peak structure 
found in BSCCO J?] may be due to the coupling 
to one optical branch with two van Hove singu- 
larities at k z = and k z = ir/d. 

The dispersion of phonons also leads to a cou- 
pling (phase-locking) of Josephson oscillations 
in different resistive junctions p|,^3[. A phase- 
locking in a stack of Josephson junctions is im- 
portant for applications of such systems as high- 
frequency mixers and detectors. 

2. Vortex motion and phonons 

In the case of long junctions and in the presence 
of an external magnetic field applied parallel to 
the layers we have to account for the variation of 
the phase 7„(x, t) along the layers (x-direction). 
The result is a generalization of the coupled sine- 
Gordon equations derived in fTTj fi~2] : 



dlln{x,t) 
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J n = sin 7„ + 
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(8) 



(9) 



The characteristic lengths for the variation of the 
phase along the layers are calculated from 
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where A c , \ a b are the penetration depths of mag- 
netic fields polarized in c-direction and parallel to 
the ab-planes, respectively. They are related to 
the corresponding plasma frequencies A c = c/loj, 
A a 6 = c/^ab, where c is the velocity of light. 
For BSCCO realistic values are A a fc = 170nm, 
d = 1.5nm, A c = 150yum, which gives a Josephson 
penetration depth of Xj ~ lum. 

In the derivation of Eq. (g) terms of the form 
7/w 2 b , which are small compared to 7 for the 
frequency range considered, have been neglected. 
Furthermore in the coupling to phonons only po- 
larization effects with polarization in c-direction 
excited by fields in c-direction (expressed by the 
dielectric function e^(k x , k z ,uj)) have been taken 



into account. It can be shown |[13| , that the ne- 
glect of other polarizations is a good approxima- 
tion for wave- vectors with k x d <C 1, which is well 
fulfilled for k x < 1/Aj. On the other hand, k z is 
not small and k z 3> k x in general. 

In the following numerical calculations we con- 
sider only local polarization effects, i.e. we ne- 
glect the dispersion of phonons. This will be suf- 
ficient to demonstrate the principle effects. 

2.1. Boundary conditions 

Assuming a stack with N Josephson barriers 
with n = 1 . . . N + 1 superconducting layers and 
two normal contacts n — 0, N + 2 this set of equa- 
tions (ph holds for n = 1 . . . N. In the barriers 
connecting the superconducting layers with the 
normal contact J n has to be replaced by the nor- 
malized bias current Jb = jb/jc- It is useful to 
incorporate this boundary condition into the set 
of equations by subtracting Jb for each term writ- 
ing J n := J n - J b then 



®lln - tj Jb 



X 2 / n 



A^( J " +1 



+ Jn-l) (10) 



which has to be solved with the boundary condi- 
tions J n = for n = 0, N + 1 for a finite stack. 

The magnetic field (in y-direction) which 
causes the variation of the phase along the lay- 
ers enters explicitly the boundary condition for 
dxln at the edges of the stack. It consists of the 
external field B ext and the field generated by the 
bias current. If we neglect the latter we have 



dx'ynix = 0) = d x j n (x = L x ) 



2ed_ 



-B„ 



17(11) 



and we may drop the small contribution Jt/A^ on 
the left side of Eq. (10). Then in the case of a 
constant phase along the layers one recovers the 
RSJ equation J n = 0. 

2.2. Analytical solution 

The set of coupled sine- Gordon equations can 
be solved numerically. Here one makes the gen- 
eral ansatz for the phases in the different layers: 

M 

7„(x,t) = T n (t)+-qx+2_\ 5"fn(m,t)cos{— — )(12) 



m— 1 



L x 



using a Fourier expansion for the spatially oscil- 
lating part of the phases and splitting off the term 



rjx which takes care of the average increase of the 
phase difference due to the applied magnetic field. 
With help of this ansatz a coupled set of differen- 
tial equations for the components 8j n (m,t) can 
be derived and solved with Runge-Kutta tech- 
niques. 

Approximate analytical solutions are possible 
by setting 

r n {t) = e n + ujt + s ln (o,t) (13) 

and linearizing the equations (pi) with respect to 
small oscillating terms S^/ n (m,t). This ansatz is 
well justified for large magnetic fields where the 
magnetic flux penetrates the stack almost homo- 
geneously. Therefore there is a voltage drop V n = 
(jn{t))fi/(2e) — u)h/(2e) over each junction. The 
oscillating part describes standing waves which 
oscillate primarily with the same basic Josephson 
frequency u>. 

Using in addition a Fourier expansion of the os- 
cillating parts in c-direction and keeping only the 
lowest harmonics in u we arrive at the following 
expression for the dc-current 

J 6 = ^-2^|-^^Im(P(fc £C ,A ;z ,c)) (14) 
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P(k x ,k z ,u) 



\I(k x ,r,)\ 2 \p(k z )\ 



£ ph 



(k x ,k z ,uj) + i 



k x c \k z 
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where the sum goes over the discrete values of 
k x = mTr/L x , k z = wr/((iV + l)d), < \n\ < N. 
The denominator contains the phonon dielec- 
tric function and the characteristic velocity 



(15) 



c(fc,)=c/yi-2^(coB(*»d)-l). 

Resonances are expected at frequencies where the 
real part of the denominator vanishes. There size 
depends on the weighting functions 



p{k z 



Z^' 



-ik z n 



(16) 



which contains a Fourier transformation of the 
static phase distribution and the function 



^^x Jo 



(17) 



depending on the magnetic field, which is peaked 

at rxix — Tj • 

3. Single contact in a magnetic field 

For a single contact, N — 1, the value of k z 
is fixed to k z — it/ (2d). Then the characteris- 
tic (Swihart-) velocity is c = Co/yl + 2\ 2 ab /d 2 = 
ujjXj, and the resonance frequencies in the I-V 
curve arc determined from 



Re 



£(«) - ( k A 2 



= o. 



(18) 



i>ii 



For a constant e^" structures appear 
current-voltage characteristic at W res ,m 



the 



m7rc/( 1 



These are the well-known Fiskc 



steps [pi. They correspond to the excitation of 
standing electromagnetic waves in the Josephson 
junction of length L x . The largest amplitude is 
obtained for wave- vectors k x ~ r\. Here the ve- 
locity of fluxons equals the phase velocity of the 
electromagnetic waves. In the case of a very long 
junction the Fiske steps merge into one flux flow 
branch which is peaked at V — fi.w ros /(2e) = 

MB^/y/^ (Eck-peak g|). 

In order to discuss the influence of phonons 
we use here for simplicity a dispersionless opti- 
cal phonon band with the dielectric function 



c ph/ 



e£M = £ o 



mi 



(19) 
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The spectrum of resonances as function of the dis- 
crete values of k x — rrnr / ' L x is shown in Fig. 1. 
One obtains two branches: for small A^-values the 
lower branch corresponds to the propagation of 
electromagnetic waves, while the upper branch is 
phonon-like. The lower branch ends at the (trans- 
verse) phonon eigenfrequency wto- The upper 
branch starts at the zero of the dielectric func- 
tion, i.e. at the longitudinal eigenfrequency wlo- 
The parameters used to calculate the dispersion 
shown in Fig. 1 are adapted to TBCCO. 

In Fig. 2 the result for the I-V curve is shown 
for three different magnetic fields. The figures 
compare numerical (Rungc-Kutta) with analyti- 
cal results. Note that the magnetic field selects 
the /c^-value where the strongest resonance oc- 
curs: fc™ ax = r\ = (2ed/fr)B cxt . In order to show 




Figure 1. Resonance frequencies for a single 
Josephson contact with one dispersionless phonon 
mode. 



the phonon effects more clearly a large separation 
between wto and wlo has been chosen. Further- 
more for numerical reasons a comparatively small 
value of the McCumber parameter f3 c — 50 has 
been used, while (3 C = 500 would be more real- 
istic (/3 C = w^/wj with uu c = 2edj c /(Ha)). For 
B cx t = only one resonance occurs at ojlo cor- 
responding to the subgap-resonance discussed in 
Sec. [|. For finite B ext one finds two groups of res- 
onances corresponding to the two branches in Fig. 
1. The fine-structure is due to Fiske-resonances 
in the stack of finite length L x . With increas- 
ing field strength the upper peak shifts to higher 
frequencies. The lower peak approaches the TO- 
frequency while loosing weight. In all cases there 
is a gap with no resonances for frequencies be- 
tween wto and ^lo- 

Generally the agreement between numerical 
and analytical calculations is good, in particu- 
lar, concerning the position of peaks. The agree- 
ment in the height of the peaks can be improved 
by going beyond the linear approximation in the 
oscillation amplitudes. The numerical calcula- 
tions show the same hysteretic behavior as the 
experimental results if the bias current density is 
changed. 
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Figure 2. I-V curves for a single contact for three 
different magnetic fields. Bo = <S>o/(dL x ) corre- 
sponds to a magnetic field with one flux quantum 
per junction. 



4. Several contacts in a magnetic field 

Numerical calculations can also be performed 
for larger stacks. Analytical calculations are no 
longer possible for N > 2 without further as- 
sumptions on the relative phases 9 n for the dif- 
ferent junctions. In Fig. 3 we show results for the 
positions of the resonances for a stack with N = 4 
junctions for all possible values of k z — nn/(bd). 
The multiple branches correspond to different val- 
ues of the characteristic velocity c which depends 
strongly on k z . 

In Fig. 4 we show results for the I-V curve for a 
large applied magnetic field. Here we find a vor- 
tex lattice which is moving due to the bias cur- 
rent. Both numerical and analytical results are 
shown as function of frequency which corresponds 
to the dc-voltage drop over a single junction. For 




Figure 3. Resonance frequencies for a stack of 4 
contacts. 



our analytical calculations of the I-V curves we 
have superimposed results calculated with all pos- 
sible values of k z . The comparison with numerical 
results in Fig. 4 demonstrates that modes with all 
possible fc z -values will be excited by cycling the 
bias current. The group of peaks at the highest 
frequency correspond to electromagnetic modes 
with the smallest k z value. The middle group 
contains phonon like excitations while the lowest 
group are electromagnetic excitations of Kleiner 
modes. Again a gap appears between wto and 
cjlo i n accordance with the dispersion curves in 
Fig. 3. 

5. Comparison with experimental results 
and conclusions 

The main results of the preceding section are: 
a) a shift of phonon resonances in a parallel mag- 
netic field, b) Fiske resonances in the frequency 
range of phonons are no longer equidistant, c) 
the flux-flow voltage is no longer proportional to 
the magnetic field and has a gap between wto 
and uJho- In order to observe these effects exper- 
imentally a strong magnetic field has to be ap- 
plied (> 3T). In most experiments |7J the applied 
field has been much lower, and no shift of phonon 
resonances has been observed. A further require- 
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Figure 4. I-V curves for a stack of 4 contacts in 
high magnetic field. 



ment is that the field is applied strictly parallel to 
the layers in order to avoid vortex pinning from 
inhomogeneities. The fact that in the cited ex- 
periments neither Fiske resonances in the fre- 
quency range of phonons nor a flux flow branch 
could be observed is an indication of vortex pin- 
ning. Recently the flux flow voltage has been 
measured for a stack of 30 junctions in BSCCO 
flq| . Deviations from the linear field dependence 
together with anomalies in the frequency range 
of an optical phonon have been observed, which 
supports the present model for the interaction be- 
tween flux flow and phonons. 
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